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A  one -dimensional,  dynamic,  thermomechanical  model,  which  includes 
nonlinear  inelastic  deformation,  internal  heat  generation  (strain-heating), 
temperature  dependent  material  properties,  thermal  expansion  and  thermoelastic 
coupling,  is  considered  for  a  uniform  thin  bar  subjected  to  mechanical  or 
thermal  disturbance.  A  nonlinear  Maxwell  material  is  examined  in  this  model 
and  special  attention  is  focused  on  the  temperature  change.  By  solving  a 
nonlinear  problem,  it  is  found  that  a  thermal  instability,  called  thermal 
runaway,  may  result  due  to  the  mutual  feedback  between  strain-heating  and  the 
temperature  dependent  inelastic  material  properties.  Neglecting  this  important 
phenomenon  may  lead  to  unexpected  material  failure./SA  linearizing  perturbation 
study  then  shows  that  the  occurence  of  this  instability  depends  on  the  choice 
of  the  material,  the  steady  state  values  of  stress  and  temperature  and  on  the 
characteristic  length  of  the  bar,  rather  than  on  the!  magnitude  or  the  form  of 
the  disturbance.  It  is  also  found  that  thermal!  expansion,  inertia  and 
therm oelastic  coupling  have  relatively  minor  effect  on  this  instability  for  the 
problem  under  consideration.  Aluminum  is  taken  as  4n  example  for  numerical 
demonstration. 
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1. 


It  is  well-known  from  thermodynamics  that  inelastic  deformations  may 
result  in  energy  dissipation  (Boley  and  Veiner  I960,  Malvern  1969).  The 
conversion  of  dissipative  mechanical  energy  into  heat  is  known  as  strain¬ 
heating,  and  this  phenomenon  has  been  studied  theoretically  by  Allen  (1985)  for 
viscoplastic  materials  and  by  Tauchert  (1967a)  for  viscoelastic  materials.  It 
has  also  been  studied  e xperimental ly  f  or  various  ma teri als  (Tauchert  1967b, 
Dillon  1962a,b,  1966)  and  significant  temperature  increases  have  been  observed, 
e.g.  a  500°K  temperature  rise  was  measured  (Dillon  1962a).  However,  few 
comprehensive  studies  have  been  carried  out  to  explore  the  conditions  under 
which  a  significant  temperature  rise  may  result.  It  is  therefore  our  objective 
here  to  include  a  wide  spectrum  of  thermal  effects,  such  as  thermal  expansion, 
thermoelastic  coupling,  temperature  dependence  of  (inelastic)  material 
properties  and  strain- heating,  in  a  one-dimensional,  dynamic,  thexmomechanical 
model.  Then  based  on  this  rather  general  model,  a  mathematical  study  of  the 
effect  of  mechanical  or  thermal  disturbance  is  given  for  a  nonlinear  Maxwell 
material  with  temperature  dependent  viscosity.  This  type  of  material  is  chosen 
because  of  the  simple  form  of  the  constitutive  law  and  because  of  its 
widespread  use  in  the  literature  to  approximate  the  behavior  of  some  polymers 
and  metallic  materials  at  elevated  temperatures. 

An  important  phenomenon  known  as  thermal  runaway,  which  is  an  unstable 
feedback  process,  is  observed.  This  thermal  instability  has  been  studied 
analytically  by  Schapery  (1964,  1965),  Huang  and  Lee  (1967)  and  experimentally 

by  Schapery  and  Cantey  (1966)  for  linear  viscoelastic  rods  subjected  to  cyclic 
loading.  It  may  occur  as  the  temperature  rise  due  to  strai  nr  heating  results  in 
a  reduction  of  the  inelastic  material  properties  (e.g.,  the  viscosity 
coefficient  in  the  case  of  Maxwell  material),  Le.  the  material  softens,  which 
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in  torn  generates  greater  inelastic  strain  and  str ai n-hea ting.  It  is  also 


found  that  thermal  expansion,  thermoelastic  coupling  and  inertia  have  a 


relatively  minor  effect  on  the  occurence  of  this  instability.  The  thermal 


runaway  phenomenon  has  been  discussed  in  the  geophysics  literature  (e.g.  Brun 


and  Cobbold  1980,  Cary  et  al.  1979  and  Van  et  al.  1986),  but  has  not  received 


much  attention  in  the  engineering  community.  This  is  probably  due  to  the  fact 


that  thermal  dissipative  mechanisms,  such  as  convection,  are  often  present  in 


conventional  engineering  applications  such  that  thermal  runaway  is  not  a 


serious  problem.  However,  in  the  absence  of  such  thermal  dissipative 


mechanisms  (e.g.,  in  space)  thermal  runaway  could  lead  to  unexpected  failure. 


In  this  paper,  we  will  study  a  thermo-mechanical  system  of  governing 


eouations  for  a  uniform  thin  bar  with  one  end  fixed  and  insulated,  and  with  the 


other  end  subject  to  a  mechanical  or  thermal  disturbance.  Thermal  expansion 


and  thermoelastic  coupling,  which  axe  well-known  effects  in  coupled 


thermoelasticity  (Boley  and  Veiner  1960)  will  be  included  in  the  material 


constitutive  law  and  in  the  energy  equation,  respectively.  Furthermore,  the 


energy  equation  will  contain  the  strain-heating  term.  These  two  equations. 


together  with  the  conservation  of  momentum  and  strain-displacement  relation. 


form  a  set  of  coupled  nonlinear  partial  differential  equations  which  govern  the 


response  of  the  system.  Particular  emphasis  will  be  placed  on  the  evaluation 


of  the  temperature  changes. 


This  nonlinear  problem  is  first  solved  in  Section  2  with  the  assumption 


that  the  load  is  applied  so  slowly  that  the  inertial  effect  can  be  ignored. 


The  resulting  problem  is  decoupled  mathematically,  and  thus  the  temperature 


field  can  be  determined  by  an  iterative  procedure  with  the  use  of  a  Green's 


function.  The  solution  to  this  nonlinear  problem  clearly  demonstrates  the 


thermal  runaway  phenomenon.  Next,  employing  the  standard  perturbation 


technique,  linearized  equations,  which  are  amenable  to  a  stability  analysis, 
are  obtained  and  studied  in  Section  3.  Solutions  are  obtained  by  the  method  of 
separation  of  variables,  and  then  the  stability  analysis  is  performed. 
Analytical  solutions  under  quasi-static  conditions  are  obtained  both  for  the 
case  of  a  mechanical  disturbance  and  for  a  thermal  disturbance.  Numerical 
results  are  presented  for  an  aluminum  bar,  and  some  concluding  remarks  are 
given  in  Section  4. 
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2.  .The 

Consider  &  thin  bar  length  L  with  one  end  (x  =  0)  fixed,  and  with  the 
other  end  (x  =  L)  subject  to  an  axial  stress  load  p(t)*,  the  bar  is  initially  at 
a  uniform  temperature  To.  The  constitutive  law  for  a  very  wide  class  of 
inelastic  materials  may  be  expressed  as 


Q( e )  =  P(o)  +  R(T)  +  g(o,a*.  T) 


where  e(x,t),  cr(x.t)  and  T(x,t)  are  respectively  the  normal  strain,  normal 
stress  and  temperature.  These  terms  in  general  vary  with  the  axial  coordinate  x 
and  the  time  t.  In  equation  (1),  P,  Q,  and  R  are  linear  differential  operators 
in  time  and  g  is  a  general  nonlinear  fucation,  which  are  chosen  to  characterize 
the  viscoelastic  or  viscoplastic  behavior  of  the  material.  For  a  fairly  large 
class  of  such  materials,  equation  (1)  can  be  simplified  to 

It  *  i  If  *  -  S ♦  (2> 


The  first  term  and  the  second  term  on  the  right  hand  side  of  equation  (2) 
represent,  respectively,  the  instantaneous  linear  elastic  response  and  the 
thermal  expansion,  which  are  characteriz ed  by  Young's  modulus,  E,  and  by  the 
linear  thermal  expansion  coefficient,  a.  The  last  term  in  equation  (2)  is  the 
inelastic  strain  rate  which  depends  directly  on  cr  and  e  with  the  dependence  on 
T  occurring  indirectly  through  the  temperature  dependence  of  the  material 
properties.  Various  forms  for  f  have  been  postulated  for  different  materials 
and  for  different  loading  situations  (e.g.  see  Cristescu  1967,  Cernocky  and 
Krempl  1980)  . 

The  conservation  of  energy  equation  consistent  with  equation  (2)  (e.g.  see 
(hang  and  Cozzarelli  1977  for  nonlinear  thermoviscoelastic  materials)  is  given 
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vhere  k,  p  and  c  are  the  thermal  conductivity,  density  and  heat  capacity  of  the 
material,  respectively.  The  second  term  in  equation  (3)  is  known  as  the 


strain—  hea ting  term,  whereas  the  third  term  represents  the  thermoelastic 
coupling  effect.  Note  that  heat  flux  through  the  lateral  surface  has  been 
neglected,  i.e.,  the  lateral  surface  of  the  bar  is  assumed  to  be  insulated. 
Also,  it  is  assumed  that  all  of  the  mechanical  energy  during  the.  deformation 
process  is  transformed  into  heat.  Equations  (2)  and  (3),  together  with  the 
conservation  of  momentum 


i 
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da 


dx 


and  the  strainrdispl  a  cement  relation 


(4) 


e 


dn 

dx 


(5) 


where  u(x,  t)  is  the  axial  displacement,  form  a  set  of  equations  which 
completely  describe  the  response  of  the  system. 

In  particular,  a  nonlinear  Hqxwell  model  with  temperature  dependent 
viscosity  will  be  considered  here.  For  this  type  of  material,  equations  (2-4) 
can  be  rewritten  with  the  use  of  equation  (5)  ss 


iz 

dx 


Ida 

Edt 


+  A(T)  on  + 


where 


(6) 


A(T)  =  Aoe~B/T 


6 


is  the  reciprocal  viscosity 


and 


2 

k  ^  +  A(T)  on+1 

3xZ 


-aT 


da 
o  dt 


da  dv 
dx  "  P  at 


(7) 


(8) 


The  variable  v(x,  t)  is  the  axial  velocity,  and  Aq,  B  and  n  are  respectively  the 
pre-exponential  constant  in  the  reciprocal  viscosity,  the  creep  activation 
constant  and  the  stress  power.  Equations  (6-8)  represent  a  system  of  coupled, 
nonlinear  partial  differential  equations  for  a  nonlinear  Maxwell  material  in 
terms  of  the  dependent  variables  v,  a,  and  T.  In  general,  it  is  difficult  to 
solve  these  equations  because  of  the  presence  of  strong  nonlinearities  and 
coupl ing. 

Now,  let  us  consider  a  special  case  in  which  the  characteristic  time,  tQ, 

e 

of  the  applied  stress  is  large  enough  so  that  inertia  effects  can  be  neglected. 

It  follows  immediately  from  equation  (8)  that  the  stress  field  is  uniformly 

distributed  and  simply  equal  to  the  applied  stress  at  the  free  end.  As  a 

result,  the  problem  is  uncoupled  and  thus  the  temperature  and  velocity  fields 

can  be  obtained  successively.  As  an  example,  it  is  assumed  that  the  bar  is 

insulated  at  the  fixed  end  and  maintained  at  a  constant  temperature,  T  ,  at  the 

o 

free  end.  The  Green's  function  associated  with  the  auxiliary  problem,  in  which 
the  known  heat  source  term  in  equation  (7),  (i.e.,  the  second  and  the  third 
terms)  are  suppressed;  is  simply 

m 

G(x, tlx ' ,t)  =  ^  £  exp [-  ~  P^(t-c)]  cosPn  x  cosPBx'  (9) 

where  8  =  xv-— 1 -  .  If  an  initial  guess  for  the  solution  is  available,  say 

m  2L> 
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(r,  t) ,  then  the  following  integrals  (e.g.  see  Ozisik  1980) 


T=t  x'=L 


T(i+1)  (x.t)=  —  f  f  G(x, tlx'.r)  [A<T(i)(x',-c))p(-c)n+1 

x=0  x  -0 


-aT  dx'dr 

O  OT 


i  =  0,1.2.. 


(10) 


can  be  carried  oat  successively  by  ntu&erical  methods  and  will  converge  to  the 
desired  solution.  The  rate  of  convergence  of  this  iterative  process  depends 
greatly  on  the  initial  guess.  Usually,  the  solution  for  the  completely 
insulated  case  (which  can  be  easily  calculated)  will  be  a  reasonable  choice  for 
the  initial  guess. 

The  solution  for  an  aluminum  bar  subjected  to  an  applied  stress  of  the 

form 

p(t)  *  <xo  (l-e-t/to)  (11) 

at  the  free  end  is  shown  in  Figs.  1  and  2.  The  material  properties  and  other 
necessary  parameters  are  given  in  Table  l;  the  creep  data  were  taken  from 
Garafalo  (1967),  Walter  and  Ponter  (1976).  Fig.  1  shows  the  temperature 
increase  with  time  at  x  =  0,  L/2  and  3L/ 4,  whil  e  Fig.  2  show s  the  temperature 
distribution  over  x  at  t  =  250,  500,  750  and  1000  seconds.  These  two  figures 
clearly  indicate  that  without  appropriate  heat  dissipating  mechanisms,  the 
temperature  can  increase  without  bound  leading  eventually  to  failure. 

It  is  important  to  note  that  although  the  above  problem  has  been  decoupled 
mathematically  by  considering  the  quasi-static  case,  it  is  still  coupled 
physically  because  the  strain-heating  feedback  into  the  temperature  dependent 
viscosity  has  already  been  imbedded  in  the  energy  equation  and  the  material 
constitutive  law.  Also,  it  is  found  in  this  example  that  thermal  expansion  and 
thexmoelastic  coupling  have  a  relatively  minor  effect  on  the  temperature  rise. 


3 •  The  Linearized  Problem  -  Stability  Analysis 

It  has  been  demonstrated  in  the  previous  section  that  the  temperatnre  can 

increase  significantly  dne  to  strain-heating  and  may  lead  to  eventual  failure. 

To  prevent  this  from  occurring,  a  cooling  process  can  be  designed  to  extract 

the  heat  generated  internally  and  thus  control  the  temperature. 

Mathematically,  this  can  be  accomplished  by  introducing  a  heat  sink  term 

=  -A(TQ)p  (t  )n+  *»  into  the  energy  equation.  As  a  result,  the  system  will 

reach  a  steady  state  in  which  T  =  T  ,  a  =  a  and  v  =  v  x/L  where  v  = 

o  o  o  o 

A(To)CTonL.  In  this  section  we  assume  that  the  system  is  initially  in  this 
steady  state  condition. 

If  the  system  is  then  subjected  to  a  small  additional  external  disturbance 
in  stress  or  temperature  at  x  =  L,  this  disturbance  will  result  in  additional 
increments  of  stress,  a,  temperature  ¥,  and  velocity,  v.  It  is  assumed  that 
the  magnitudes  in  the  steady  state  are  much  larger  than  the  magnitudes  of  the 
increments  in  the  perturbed  state,  and  we  may  then  write 


T  =  T  +  T,  T  >>  T 
o  o 


a  -  aQ  +  a,  oq  >>  CT 


V  =  Vo  +  v*  vo  »  v 


(12) 


Substituting  equation  (12)  into  equations  (5-8),  subtracting  the  relation  for 
the  steady  state  (with  sink  Sq  added  to  equation  (7)),  retaining  the  linear 
terms  in  the  increments  (as  in  Cozzarelli  et  al.  1970),  and  then  introducing 
the  nondimensional  quantities 


we  obtain  the  following  linearized  equations: 


u 


3v  1  da  _  3T 

3z  =  En  +  D0+7lT+a3T 


da 


3  v 


3x  p  3t 


k  i!|  +  (n+l)  a  +  qT  -  a  ff  =  pc  ff 
3x 


(13) 


(14) 


(15) 


Note  that  the  overbars  (  )  have  been  deleted  in  equations  (13-15)  for  the  sake 
of  convenience  of  notation.  The  associated  boundary  and  initial  conditions  are 
as  follows: 


v  =  0.^  =  0 
'  3x 


at  x  =  0 


«a«1,  T  =  Tj  at  x  =  1 
T  =  a  =  v  =  0  at  t  =  0 


(16) 


where  o ^  =  0  for  a  thermal  disturbance  and  =  0  for  a  mechanical  disturbance. 

It  can  be  seen  that  the  method  of  separation  of  variables  may  be  employed 
if  we  separate  equations  (13-16)  into  two  problems  (see  Street  1973),  i.e.  a 
steady-state  problem  with  inhomogeneous  boundary  conditions  and  a  transient 
problem  with  homogeneous  boundary  conditions.  Solving  these  two  problems  and 
superposing  solutions,  we  obtain 


T  =  Ti  +  1  Pa(t)  cosPmx 
bf=1 


(17) 
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c  "  ®i  +  1  9m<t)  cosPDx 
m=l 


v  =  vi  +  1  *T»(t)  sin^j 
m=l 


where 


( 2a- 1 )  tt 

Pm  =  2 


The  temporal  functions,  Pm(t),  qn(t)  and  rm(t)  are  to  be  determined.  In 


equations  (17-19),  T^,  <y^  and  are  the  following  solutions  for  the  steady 


state  problem: 


~  =  a+1  ,  cos  fee 
i  q  cos£  *  C1 


®i  =  °1 


▼ i  =  «1  t-x+(tt+l)  ] 


§cos§ 


for  a  mechanical  disturbance,  and 


_  _  T  cog& 
Ai  "  ll  co s 5 


Oi  =  ° 


sintx 
vi  =  Tln  i|cos$ 


for  a  thermal  disturbance,  where 


« -  qr 


1  * 1 


«  1  *  *  t  '  ♦  1  *  *  n 


$5 
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Substituting  equations  (17-19)  into  equations  (13-15).  multiplying  by  cos^x, 

sinS  x  and  cosB  x  (  =  1,2,....),  respectively,  integrating  tbe  resulting 

Xj  x* 

equations  over  the  space  domain  (x  B  0  to  x  =  1)  and  invoking  the  orthogonality 
properties,  the  quantities  pm(t),  qffl(t),  and  rm(t)  equations  (17-19)  can  be 
solved  approximately  for  m  =  1,2,..,M.  Accordingly,  we  obtain  the  formulation 

Z  -  A  Z  (20) 
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where  I  is  the  identity  matrix,  and  C  and  D  are  the  diagonal  matrices 


C  =  diag(rj-k3az)  m  =  1,2,..,M 


D  =  diag({Jffl)  1,2,..,] 

As  initial  conditions,  we  have 


t>(„,  . 


,.(0)  =  (-l)"*1  2r  „ 

B 


r„(0)  -  (-1)"*1  2(^_+ 

B  2  S  -B  2 


for  the  mechanical  disturbance,  and 


,  nrf  1  .  „  2  1 


(-1)  2B  r  T. 

1 


q  (0) 


r  (0)  =  (-1) 
SI 


nrf  1  _ 2t 


52-Pb2 


for  the  thermal  disturbance.  To  obtain  accurate  results,  a  reasonable  number 
of  terms  in  equations  (17-19)  (e.g.  M2_6)  should  be  considered.  However,  this 
will  introduce  difficulties  in  the  numerical  integrations,  because  matrix  A  is 
of  relative  large  order.  Also,  because  of  the  inertia  effect  the  matrix  A  is 
stiff,  and  thus  a  very  small  time  step  has  to  be  taken  for  each  integration. 
These  two  facts  make  it  very  difficult  to  obtain  long  time  solutions  by  direct 
integration.  Alternatively,  transform  methods  (e.g.,  see  Meirovitch  1980,  Chen 
1984)  may  be  employed  to  solve  equation  (20). 

If  as  in  the  nonlinear  problem,  we  again  neglect  the  inertia  effect,  the 
solution  of  equations  (13-15)  may  be  found  with  the  use  of  variable 
transf or® a ti one  (Ozisik  1980,  Carslaw  and  Jaegger  1959)  as 

T(x,  t)  *  -■ n— 1-  o.  5  I  (l-e***)  cosB  x  (21) 


for  the  mechanical  disturbance  and 

T<*-*>  *  Hi?  ti  *  7177 eV 

«^l  5  -Pm 

for  thenaal  disturbance.  In  equations  (21)  and  (22), 


*«v«* 


4'  f  f ,  v  *  l' 


X  =  —  (q-k£  2) 
m  pc  1  MB 

It  is  interesting  to  note  that  the  same  XB  appears  in  equations  (21)  and  (22). 
This  indicates  that  the  perturbations  in  temperature  will  follow  the  same 
general  pattern,  regardless  of  the  forms  and  the  magnitudes  of  the 
disturbances.  Furthermore,  for  this  quasi-static  solution  to  be  stable,  it  is 
required  that 


(23) 


or  in  dimensional  fon 


BL2A  e_B/Too  ttfl 
o  _ o 

To2k 


(24) 


Using  Jhe  values  of  the  parameters  listed  in  Table  1,  it  is  found  that 
conditions  (23)  or  (24)  axe  violated.  Thus,  the  quasi-static  system  with  these 
values  for  the  parameters  is  unstable. 

Let  us  now  return  to  a  consideration  of  stability  for  the  dynamic  problem. 
Substituting  equations  (17-19)  into  equations  (13-15),  and  combining  the 
resulting  equations  into  one  equation  in  terms  of  p  gives 


•lPm  +  *2 Pm  +  *3Pm  +  a4Pa 


where 


(25) 


al  ■  e  <PC“«2e) 
kin 

*2  “  E  ^m  ~  E  +  pcn  +  a^n+D  “  a1\ 
*3  “  kPa2n  +  cPB2  +  H 

%  -  i  <kP»2  -  p» 


(26) 
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lie  overdots,  (*),  represent  derivatives  with  respect  to  the  time.  By  applying 
the  Routh-Hurw  i  tz  criterion  (e.g.,  see  Kuo  1982),  one  root  (i.e.  m  =  1)  with  a 
positive  real  part  is  found  for  aluminum  using  the  parameters  listed  in  Table 
1.  This  observation  indicates  that  inertia,  thermal  expansion  and 
thermoelastic  coupling,  which  were  not  considered  in  obtaining  equations  (21) 
and  (22),  do  not  contribute  to  the  instability  found  earlier.  Thus,  it  may  be 
concluded  that  the  mechanisms  which  may  cause  thermal  instability  (runaway)  are 
the  temperature  dependence  of  inelastic  material  properties  and  the  strain¬ 
heating.  This  conclusion  can  also  be  obtained  by  letting  B  =  0  (or  q  *  0  in 
nondimensional  form)  in  equation  (26).  The  resulting  coefficients  in  equation 
(25)  are  then  all  positive  (see  Chang  and  Cozzarelli  1977  for  a^  >  0)  and 

2 

*2a3  “  ala4  =  <£  n  +  f”*  kpB4  +  t<>cn2k  +  P®2  +  o(n+1)  ** 

+  o(nt-l)c]p„2  >  0 

n 

Using  the  Routh-Hurw  i tz  criterion  again,  the  system  is  seen  to  be  stable  in 
this  case.  Ihus,  if  the  temperature  dependence  of  the  material  properties  is 
not  considered,  the  thermal  runaway  phenomenon  will  not  be  observed. 

It  is  important  to  note  that  in  accordance  with  condition  (23)  the 
occurrence  of  thermal  instability  depends  on  one  nondimensional  parameter,  g, 
whose  value  is  determined  by  the  material  properties  AQ*  B,  n,  k,  the  steady 

i 

state  values  T  a Q,  and  the  length  of  the  bar  L.  Thus,  in  the  linearized 
problem  the  occurrence  of  runaway  is  not  affected  by  the  magnitude  or  form  of 
the  disturbance.  The  stability  relationship  between  Tq  and  <sQ  (Le.  condition 
(24))  for  an  aluminum  bar  of  length  1  meter  is  shown  in  Fig.  3.  Also  shown  in 
Figs.  4  and  5  are  the  temperature  increments  due  to  mechanical  and  thermal 
disturbances,  respectively,  where  Oj/a0  and  Tj/Tq  »re  given  at  the  end  of  Table 


4 .  Cone lad ina 


A  one-dimensional,  dynamic,  thermomechanical  model,  which  includes 
nonlinear  inelastic  deformation,  internal  heat  generation  (strain-heating), 
thermal  expansion,  thermoelastic  coupling  and  temperature  dependent  material 
property,  is  considered.  A  detailed  study  is  given  for  a  special  material 
model  (nonlinear  Maxwell),  which  approximately  characterizes  some  polymers  and 
metallic  materials  at  elevated  temperatures.  For  a  quasi-static  nonlinear 
problem,  significant  temperature  increases  are  found  for  aluminum  by  i^ans  of 
an  iterative  numerical  solution.  Further  investigations  by  a  linearized 
perturbation  technique  shows  that  a  thermal  instability  (runaway)  may  occur, 
depending  on  the  choice  of  material,  steady  state  values  of  stress  and 
temperature  and  the  length  of  the  bar,  but  not  on  the  magnitude  or  the  form  of 
the  disturbance.  The  phenomenon  of  thermal  runaway  is  the  result  of  mutual 
feedback  between  strain- heating  and  the  temperature  dependent  (inelastic) 
material  properties.  If  either  effect  is  not  included  in  an  analysis,  the 
possibility  of  thermal  runaway  will  not  be  explored  and  an  unexpected  failure 
may  occur.  The  effect  on  the  thermal  runaway  due  to  inertia,  thermal  expansion 
and  thermoelastic  coupling  has  been  shown  to  be  relatively  minor  for  the 
considered  one-dimensional  problem  However,  these  effects  could  be  important 
for  more  complicated  loading,  geometric  situations  and  flexible  structures 
(e.g.  trusses,  beams,  plates).  Such  cases  require  further  study. 
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Table  1  Material  Properties  for  Aluminum  and  Other  Parameters 


density  p  (Kg/m3) 

2707 

creep  activation  constant  B  (°K) 

13603 

heat  capacity  c  (J/Kg°K) 

896 

length  L  (m) 

H 

thermal  conductivity  k  (W/m°K) 

222 

steady  state  stress  a (Pa)  3 

.  5xl07 

Youngs'  Modulus  E  (Pa) 

5.52xl010 

reference  temperature  Tq  (°K) 

533 

stress  power  parameter  n 

4.55 

characteristic  time  t  (second) 

o 

10 

pre-exponential  constant  Aq 

7.836xl0"27 

mechanical  disturbance  a^/o^ 

.05 

thermal  expan.  coeff.  o  (1/°K) 

2 . 3xlo”5 

thermal  disturbance  T,/T 

1  o 

.01 

Figure 

Fig.  1 

Fig.  2 

Fig.  3 

Fig.  4 


Caption? 

Nonl  inear  probl  ea  -  temperature  increase  with  tat  x  =  0 ,  L/2.  3L/4 
for  alnainum 

Nonlinear  problem  -  temperature  variation  with  x  for  t  =  250,  500, 
750,  1000  sec.  for  aluminum 

Linearized  problem  -  stability  diagram  of  steady  state  stress  vs. 
steady  state  temperature  for  an  aluminum  bar  1  meter  long 

Linearized  problem  -  temperature  increment  due  to  mechanical 
di sturbance  - 


Fig.  5  Linearized  problem  -  temperature  increment  due  to  thermal  disturbance 


Pig.  1  Nonlinear  problem  -  temperature  increase  with  t  at  x  =  0,  L/2,  3L/4 


for  aluminum 

Nonlinear  problem  -  temperature  variation  with  z  for  t  =  250.  500. 
750,  1000  sec.  for  aluminum 

Linearized  problem  -  stability  diagram  of  steady  state  stress  vs. 
steady  state  temperature  for  an  aluminum  bar  1  meter  long 

Linearized  problem  -  temperature  increment  due  to  mechanical 
di sturbance 

Linearized  problem  -  temperature  increment  due  to  thermal  disturbance 
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Fig.  2 


Fig.  3 


Fig.  4 


Fig.  5 
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